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ABSTRACT 

We examine the luminosity function (LF) of [Oil] emission-line galaxies in the high- 
resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy dis¬ 
tribution of each galaxy, we select a sub-sample of star-forming galaxies at 0.06 < z ^ 3.0 
using the [Oil] emission line luminosity L([OII]). We confirm that the specific star forma¬ 
tion rate matches that in the GAMA survey. We show that the [Oil] LF at z — 1.0 from 
the MB II shows a good agreement with the LFs from several surveys below L([OII])=10 43 0 
ergs -1 while the low redshifts (z < 0.3) show an excess in the prediction of bright [Oil] 
galaxies, but still displaying a good match with observations below L([OII])=10 416 erg s -1 . 
Based on the validity in reproducing the properties of [Oil] galaxies at low redshift (z < 1), 
we forecast the evolution of the [Oil] LF at high redshift (z ^ 3), which can be tested by up¬ 
coming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends 
range from -3 to -2 showing minima at z = 2. The slope of the bright end evolves approxi¬ 
mately as (z + l) -1 at 2 ^ 2 while the faint end evolves as ~ 3 (z + l) -1 at 0.6 ^ z ^ 2. 
In addition, a similar analysis is applied for the evolution of [0111] LFs, which is to be ex¬ 
plored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation 
function of [Oil] and [OIII] emitting galaxies shows a rapid evolution from z = 2 to 1. 

Key words: cosmology: theory - galaxies: evolution - galaxies: luminosity function, mass 
function - stars: formation - methods: numerical. 


1 INTRODUCTION 

Understanding the star formation history of the Universe is essen¬ 
tial in decoding the evolution and mass assembly of galaxies ( Lilly 
et al.||1996[ |Madau et al.||1996[ |Hopkins||2004[ |Hopkins & Bea- 

com 2006). Observations of star-forming galaxies at high redshift 
reveal that the cosmic star formation rate (SFR) density has been 
decreasing since reaching the peak at z ~ 2 (|Madau et al.|1996| 

|Hopkins & Beacom|2006[ |Kennicutt & Evans |2012| ). Thus, inter¬ 

preting the evolution of the global SFR around this redshift (z ~ 2) 
has been one of the main interests in the astrophysics community. 
In particular, the role of self-regulated processes of star formation 
and the effect of feedback by accretion on to supermassive black 
holes (SMBH) at the centres of galaxies have been the main clues 
in deciphering the observed star formation history. 


* E-mail: kwangho.park@physics.gatech.edu (KP); 

tiziana@phys.cmu.edu (TDM); shirleyh@andrew.cmu.edu (SH) 


The SFR of an individual galaxy is estimated using several in¬ 
dicators such as nebular emission lines or the UV continuum. Both 
indicators are directly linked to the ionizing flux from young mas¬ 
sive stars, whereas the flux in the mid-infrared (MIR) or far-infrared 
(FIR) are used to trace the re-radiated emission by dust. Among 
the star formation indicators, Ha is regarded as one of the best di¬ 
rect indicators of the current SFR since the line strength of Ha is a 
good tracer of the photo-ionized gas around massive young stars (< 
20 Myr) ( |Kennicutt|1998) . However, due to the rest-frame wave¬ 
length (Ah « - 6563 A), it is challenging to use Ha for galaxies at 
high redshift (z > 1) since the observed wavelength is redshifted 
to a wavelength longer than NIR (A 0 bs > 1 pm), which is diffi¬ 
cult to observe from the ground. In contrast, [O il] AA3726, 3729 
emission line doublet can still be observed from the ground in the 
wavelength A 0 bs < 1 pm for galaxies at high redshift (z < 1.6) 


([Gallego et al.|2002||Ly et al.|2007||Gilbank et al.|2010| Zhu, Mous-| 
|takas & Blanton|2009||Ciardullo et al.|2013||Comparat et al.|2015| ). 
Since the [O II] line luminosity L([OII]) is an indirect indicator (un- 


© 2015 RAS 

















2 K. Park, T. Di Matteo, S. Ho, R. Croft, S. Wilkins, Y. Feng and N. Khandai 


like Hck), it requires an extensive calibration based on the properties 
of the galaxies such as the metallicity and stellar mass (M*) (|Kew-| 


ley, Geller & Jansen|2004||Moustakas, Kennicutt & Tremonti|2006| 


Gilbank et a l. 2010 ). At even higher redshift (z > 2), one can use 
the UV continuum, which also directly traces the ionizing photons, 
but it measures the SFR over longer timescales (~ 100 Myr). 

Over the past few decades, systematic searches for emission¬ 
line galaxies (ELGs) have been made. |Gilbank et al.| < |2010| ) ex¬ 
plored the SFR at z ~ 0.1 using [O II], Ha and u-band lumi¬ 
nosities from the deep 275 deg 2 Stripe 82 field in the Sloan Dig¬ 
ital Sky Survey (SDSS) coupled with UY data from the Galaxy 
Evolution EXplorer (GALEX) satellite. Also, the pilot survey of 
the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX; 
|Adams et al.|2011[|Ciardullo et al.|2013] > observed 284 [O II] emit¬ 
ting galaxies at z < 0.56 in 169 arcmin 2 . The HETDEX pilot 
survey is a blind integral-field spectroscopic study of four data- 
rich areas of sky: COSMOS ([Scoville et al.|2007), GOODS-N (Gr] 
|avalisco et al .12004]), MU NICS-S2 ( |Drory et al.|2001) , and XMM- 
LSS ( [Pierre et al.||2004). The Galaxy And Mass Assembly sur¬ 
vey (GAMA pBauer et af]|2013| ) also targets ~ 73, 000 galaxies 
at 0.05 < z < 0.32 with M* < 10 10 M©. The GAMA survey 
calculates the SFR using the Ha emission line and explores how 
the specific SFR (sSFR) depends on the stellar mass of the galax¬ 
ies. Recently, Comparat et al. (2015) explored the evolution of the 
luminosity function (LF) of the [O II] emitters in the redshift range 
0.1 < z < 1.65 based on the medium-resolution flux-calibrated 
spectra of ELGs with the VLT/FORS2 instrument and the SDSS- 
III/BOSS spectrograph. In the forthcoming years, there are ambi¬ 
tious surveys targeting ELGs at high redshift. The Dark Energy 
Spectroscopic Instrument (DESI) is a multi-fiber spectrograph that 
will target 18 million ELGs ( [Schlegel et al. 120091 |Levi et al.|2013) . 
DESI aims to probe the effect of the dark energy on the expansion 
history of the Universe using baryon acoustic oscillations (BAO) 
with a sky coverage of 14,000-18,000deg 2 . Another survey is the 
space-based Wide-Field InfraRed Survey Telescope-Astrophysics 
Focused Telescope Assets (WFIRST-AFTA) which will perform a 
deep infrared survey of 2400 deg 2 and obtain spectroscopic red- 
shifts of galaxies targeting 20 million Ha emitters at 1 < z < 2 
and 2 million [O III] ELGs at 2 < z < 3 ( [Spergel et al.|2013j ). 

In cosmological simulations, it is challenging to study stellar 
mass assembly of individual galaxies as it requires a high resolution 
and a large volume at the same time. A sufficiently large volume is 
needed for statistically meaningful results whereas high resolution 
properly manifests baryonic physics related to star formation and 
black hole (BH) feedback. It is also difficult to evolve a simulated 
universe up to the present (z ~ 0). Considering all these numer¬ 
ical challenges, the MassiveBlack-II simulation (MBII; |Khandai| 
|et al. 12015), which is a successor to the MassiveBlack simulation 
( |Di Matteo et al.|2012| ), is one of the unique simulations in that it 
is a state-of-the-art high resolution cosmological hydrodynamical 
simulation with a large comoving volume (c.f., |Vogelsberger et al.| 
|2014||Schaye et al.|2015| ). 


In this paper, we investigate the properties of the star¬ 
forming galaxies in the MBII simulation, which are selected by the 
[O II] and [O III] emission-line luminosity in the generated spectral 
energy distributions (SEDs). The main goal of this study is to val¬ 
idate whether the star-forming galaxies in the high-resolution cos¬ 
mological simulation can reproduce results from recent surveys for 
the [O II] emitting galaxies and make a prediction for the upcoming 
surveys that target the high redshift [O II] emitters. In Section[2]we 
explain how we select the [O II] emitting galaxies from the syn¬ 
thesised SEDs and compare the properties of the selected samples 


with observations. In Section[3] we compare the LFs from the sim¬ 
ulation and observation at 0.1 ^ z ^ 1.0, and make a prediction 
for the LF of the [O II] emitting (and [O III]) galaxies at high red¬ 
shift. We also present the evolution of the auto-correlation function 
of the [O II] ELGs at z ^ 4. Finally, we summarise and discuss the 
results in Section 0] 


2 METHODS 

2.1 The MassiveBlack-II simulation 


We use the state-of-the-art high resolution hydrodynamical simula¬ 
tion MBII ( [Khandai et al.|2015| ) performed with the hybrid version 
TreePM-SPH code P-GADGET designed for running on PETAFLOP- 
scale supercomputer facilities. The MBII contains N part = 2 x 
1792 3 dark matter and gas particles in a comoving volume of 
(100 h~ x Mpc) 3 to satisfy the needs of large volume and high res¬ 
olution at the same time. In the MBII simulation, the ACDM uni¬ 
verse is evolved from z — 159 to z — 0.06 with a high mass reso¬ 
lution of 1.1 x 10 7 h~ 1 M © for dark matter and 2.2 x 10 6 /i _1 M 0 
for gas particles with a smoothing length of 1.85 h~ 1 kpc. The cos¬ 
mological parameters used in the MBII are: erg = 0.816, spectral 
index n s — 0.968, matter density fraction of the critical density 
D m = 0.275, vacuum energy density fraction Da = 0.725, baryon 
density fraction Db = 0.046, and Hubble constant h = 0.702 from 
the 7-year Wilkinson Microwave Anisotropy Probe ([Komatsu et al.| 

[20TT ). 


The MBII includes a sub-grid model for star-forming multi¬ 
phase gas ( [Springel & Hernquist[2003| . In this model, a thermal in¬ 
stability takes place at a critical density threshold creating a multi¬ 
phase medium consisting of cold clouds in pressure equilibrium 
with surrounding hot gas. A star formation prescription is given 
by the Kennicutt-Schmidt law ( |Kennicutt|1989| >, where the SFR is 
proportional to the density of cold clouds (i.e., psfr oc p^ as where 
N = 1.5 is adopted). Gas particles are converted to star particles 
according to the star formation prescription. Star formation is then 
regulated by the supernovae feedback, which heats the surrounding 
gas and creates a self-regulated cycle. 

In the MBII, a BH is introduced as a collisionless sink parti¬ 
cle in a newly collapsing halo identified by the friends-of-friends 
halo-finder on-the-fly at a regular time interval. A seed BH with 
mass M S eed = 5 x 10 5 /i _1 M 0 is inserted into a halo with mass 
Mhaio ) 5 x 10 10 h~ 1 M Q . After seeded, BHs are assumed to 
grow at the Bondi-Hoyle rate ( [Bondi & Hoyle|1944||Bondi|1952| ) 
Mbh = 47tG 2 PooM| h (cs,oo + Vbii)- 6/2 , where v B h is the ve¬ 
locity of the BH relative to the surrounding gas, and and c S;00 
are the density and sound speed of the gas in the multi-phase state 
( Pelupessy, Di Matteo & Ciardi 20 071). The accretion rate is lim¬ 


ited by 2 x MEdd (Begelman, Yolon teri & Rees 


Rees 


2006 


Yolonteri & 


2006 where MEdd is the Eddington accretion rate. The accre¬ 


tion rate is converted to a bolometric luminosity Lboi = ^Mbhc , 
where p is the radiative efficiency and we adopt the standard value 
of 0.1 for a thin disk model ( [Shakura & Sunyaev|1973| ). It is also as¬ 
sumed that 5 per cent of the BH luminosity thermally couples with 
the surrounding gas, isotropically depositing the radiation energy 
to the gas particles within the BH kernel (64 nearest neighbours) to 
match the observed Mbb-ct relation ([Di Matteo, Springel & Hern-J 
|quist||2005| |Springel, Di Matteo & Hernquist||2005| ). The current 
BH growth model has been adopted in extensive studies (Li et al. 


2007; Sijacki et al. 2007; Di Matteo et al. 2008; Colberg & Di Mat- 


teo||2008| Croft et al.||2009[ Degraf, Di Matteo & Springel 2010[ 
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Figure 1. Examples of the spectral energy distribution of 20 star-forming 
galaxies in the range 10 < SFR (Mq yr -1 ) < 100 at 2; = 0.1 in the 
MBII simulation. The stellar continuum and nebular emission lines are gen¬ 
erated from the star formation history of each galaxy using the stellar pop¬ 
ulation synthesis code pegase.2. Note that the nebular emission lines such 
as Lya, [O 11 ] AA3726, 3729, [O ill] A5007, and Ha are distinct. 


201 l||Degraf et al.|201 l||Di Matteo et al.|20l2l|DeGraf et al.|2012| 

Chatteijee et al. 2012, Fe ng et al.|2014) . BHs are assumed to merge 
when one BH enters within the kernel of another BH with a relative 
velocity below the local gas sound speed. 


2.2 [Oil] emission-line galaxies in the MBII 

The high resolution of the MBII enables us to generate a SED of 
an individual galaxy. We adopt the previous work by |Wilkins et al.| 
< |2013b]) to generate the SEDs and we briefly explain the process in 
this section. 


2.2.7 Generating SEDs with emission lines 


The stellar population synthesis code PEGASE.2 ( |Fioc & Rocca-| 
|Volm erange 1997 1999) is used to generate SEDs of individual 
star particles as a function of the stellar mass, age, and metallic- 
ity and assuming the Salpeter initial mass function. We then cal¬ 
culate hydrogen line fluxes first using the PEGASE.2 while non¬ 
hydrogen line fluxes are estimated from hydrogen lines using 
the metallicity-dependent conversion by [Anders & Fritze-v. Al-| 
|vensleben| ( |2003| ). These nebular emission lines, reprocessed ioniz¬ 
ing radiation by interstellar medium, include Lya, [O II] AA3726, 
3729, [O III] A5007, and Ha. The rest-frame SEDs of galaxies are 
generated by integrating all the SEDs of star partic les and e missio n 
lines with a wavelength resolution of 20 A (Wilkins et al. 2013b|ai. 
For a galaxy with a stellar mass M* = 10^ h~ l M©, for example, 
^450 SEDs from individual star particles are integrated to produce 
the SED of the galaxy. Fig. [l] shows 20 randomly selected exam¬ 
ples of synthesised SEDs of star-forming galaxies with 10 < SFR 
(M© yr -1 ) < 100 at z — 0.1. Note that the nebular emission lines, 
such as Lya, [O n], [O III], and Ha, are clearly visible for these 
star-forming galaxies. 
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Figure 2.2D histogram of the [O n] luminosity and the SFR in the MBII for 
different redshifts (z =0.06, 0.1, 0.3, 1.0, 2.0 and 4.0). The dashed line in 
each panel shows the empirical relationship from |Kewley , Geller & Jansenj 
{20041, which shows a good agreement with the MBII simulation. The 
galaxies with higher SFR show a tight match while the scatter increases 
as the L([0 ii]) decreases. 


2.2.2 Selecting [Oil] emission-line galaxies from SEDs 

We select a sample of [O II] ELGs using the generated SED for 
each galaxy. The flux of [O II] emission line for an individual 
galaxy is calculated from the SED by mimicking observation. We 
apply a simple criterion LA373o/L C ont > 1 , where L CO nt is the esti¬ 
mated flux of the continuum at A3730 by averaging the fluxes at the 
neighbouring wavelengths. This method selects galaxies with posi¬ 
tive flux at the [OII] emission line wavelength. Then the luminosity 
of the [O II] emission line component is obtained by subtracting the 
estimated continuum as L ([0 II]) = La 373 o — L CO nt. We check that 
this simple method recovers the original line [O II] luminosity well 
(see Appendix |A|). 


2.2.3 Comparison with L([OII])-SFR empirical relation 

Fig. HI shows the relation between L([0 II]) and SFR of the 
[O II] ELGs at redshifts z — 0.06, 0.1, 0.3, 1.0, 2.0, and 4.0. The 
MBII simulation displays a good agreement with the empirical re¬ 
lationship by |Kewley, Geller & Jansen| ( |2004| ), 

SFR (M 0 yr' 1 ) = 1 ^ (erg s" 1 ), (1) 

which is shown as a dashed line in each panel. In general, the 
galaxies with higher SFRs show a better agreement with the em- 
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Figure 3. 2D histogram for the distribution of SFR and stellar mass (M*) 
of galaxies in the MBII at the redshift z =0.06, 0.1, 0.3, 1.0, 2.0 and 4.0, 
respectively. Most of the star-forming galaxies with M* > 10 9 r'Mg at 
z = 0.06 are not selected by the L([0 n]) while most of the star-forming 
galaxies with the same stellar mass at z = 4.0 are selected. The dashed line 
for each redshift is from |Speagle et al.| < |2014) . 



log(L[OII]) > 41.5 (z=0.06, 0.1, 0.3) 
log(L[OII]) > 41.0 (z=0.06, 0.1, 0.3) 
All SF galaxies (z=0.065, 0.1, 0.3) 

• GAMA (0.05 <z <0.32) 




log(M, h~ l [M 0 ]) 


Figure 4. Specific star formation rate (sSFR) as a function of the stellar 
mass (M*/M©) at z = 0.06, 0.1, and 0.3. Different symbols are used for 
different luminosity cuts: triangles are all star-forming galaxies, squares are 
galaxies with L([0 n]) ^ 10 410 erg s _1 , and pentagons are galaxies with 
L([0 II]) > 10 41 - 5 erg s 1 . The shaded region shows the mean sSFR with 
scatter from the GAMA survey (Bauer et al.|2013) at 0.05 < z < 0.32 as 
a function of the stellar mass in the range 9.0 < log (M* [M©]) < 11.0. 


log SFR= (0.84 =b 0.02 - 0.026 zb 0.003 xf)logM* 

- (6.51 =b 0.24 - 0.11 zb 0.03 x t) 

( 2 ) 

where t is the age of the universe in Gyr, is overplotted in Fig. [3] 
Since it is known that the MBII produces too many galaxies at low 
and high stellar mass (Khandai et al.|2015) , it is expected that the 
overabundance in stellar mass function inevitably propagates into 
the [O II] LF. The excess is obviously seen in Fig.[3]at low redshifts. 


pirical relation while the scatter increases as L([0 II]) decreases. 
Due to the mass resolution of star particles in the MBII (i.e., 
ra* as 2.2 x 10 6 /i -1 M©), shot noise starts to dominate for 
the sample of galaxies with L([0 II]) ^ 10 40 ' 6 ergs -1 (see Ap¬ 
pendix [A] for details). The current SFR is not well represented in 
L([0 II]) below this luminosity cut, which sets the lower limit of 
L([0 II]) for our current study. This luminosity cut can be trans¬ 
lated to a SFR of ~ 0.3 M© yr -1 using equation {!]>. 

Fig. [3] shows the distribution of the stellar mass (M*) and 
SFRs of all star-forming galaxies in the MBII for redshifts z =0.06, 
0.1, 0.3, 1.0, 2.0 and 4.0. In general, the SFR is proportional to the 
stellar mass M* for each redshift, and the SFR for a given stel¬ 
lar mass increases as we go to high redshift. Note that the num¬ 
ber of stellar particles is ~450 for a galaxy with stellar mass 
M* ~ 10 9 /i -1 M©. When the luminosity cut mentioned above is 
applied (i.e., SFR of ~ 0.3 M© yr -1 ), it removes a large fraction of 
galaxies at low stellar mass while the sample of star-forming galax¬ 
ies is more complete for massive galaxies. The average stellar mass 
of galaxies selected this way is approximately M* ^ 10 10 
at z = 0.06. The completeness at a given stellar mass also depends 
on redshift. For example, most of the star-forming galaxies with 
M* ~ 10 9 h~ 1 M© at z — 0.06 are not selected for our luminosity 
cut while most of the star-forming galaxies with the same stellar 
mass are selected at z — 4.0. 

For comparison, the observed relation between the SFR and 
the galaxy stellar mass from |Speagle et al.| ( |2014| ) 


2.2.4 Specific star formation rate 

In this section, we compare the SFR per stellar mass (i.e., 
sSFR=SFR AM*) of the [O II] emitters with the GAMA survey 
(Bauer et ah]|2013| ). Fig. El shows the average sSFR as a func¬ 
tion of the stellar mass in the range 9.0 < log(M* [M©]) < 
11.0 for redshifts z = 0.06, 0.1, and 0.3 with different L([0 II]) 
cuts. Triangles are all star-forming galaxies, squares are galaxies 
with L([0 II]) ^ 10 41 ° ergs 1 , and pentagons are galaxies with 
L([0 II]) ^ 10 41 ' 5 ergs -1 . The shaded region in Fig. [ 4 ] shows 
the observed average sSFR with scatter from the GAMA survey. 
The mean sSFR increases with higher luminosity cuts and the 
sSFR with the luminosity cut of 10 41 ' 5 erg s -1 matches well with 
the GAMA survey. However, this luminosity cut removes galaxies 
with M* ^ 10 10 M©. The simulation converges to the result of 
GAMA survey with increasing luminosity cut, but the halos with 
stellar mass M* < 10 10 M© are lost with this luminosity cut 
(L([Oll])^ 10 41 ' 5 erg s -1 ). 

2.2.5 Possible contamination from active galactic nuclei 

In our scheme, the SED of a galaxy does not include the emission 
lines from the hot ionized gas produced by active galactic nuclei 
(AGN). We assume that the whole [O II] emission line flux origi¬ 
nates from SFR. AGN contamination to [O II] line certainly affects 
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• Comparat et al. (2015) z=0.17 
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z=1.0 
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• Comparat et al. (2015) z ~ 0.95 
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Figure 5. Comparison of the LFs of the [O n] emitting galaxies in the MBII simulation at z = 0.1 (top), 0.3 (middle), and z = 1.0 (bottom) with the 
observations. Big circles show the [O ii]-selected sample in the MBII simulation and the small circles along with the Poisson errors show the samples 
with AGNs excluded using the method mentioned in Section |2.2.5| to avoid the possible contamination of AGNs. There is no severe contamination to the 
[O II] emission line. Red (observed) and blue (dust-corrected) solid lines at z = 0.1 and 0.3 are from HETDEX pilot survey ( [Ciardullo et al.|2013) while the 
red dashed line (observed) at z = 0.1 shows the observation from SDSS ( [Gilbank et al.|2010] >. The LFs of the MBII show an excess in the prediction_of bright 
[O II] emitters at L([Q II]) > IQ 41 - 6 erg s -1 at low redshift (z < 0.3). At z = 1.0, the LF matches well with the DEEP2 observation i Zhu, M oustakas & 
|Blanton|2009) and jComparat et al.|j2015| below L([0 ii])=10 43 0 ergs -1 . 


the estimation of the SFR of a galaxy, but observations show that 
the level of actual contamination is not high (e.g., |Bauer et al.|20~l3| 
|Zhu, Moustakas & Blanton|2009| >. 

However, we exclude galaxies that might potentially contam¬ 
inate our samples when the AGN-associated [O II] emission line 
is added to the current SEDs. We only select galaxies where the 
[O II] luminosity from star formation is larger than that from the 
AGN as L([0 ii])agn < L([0 ii])sf. From the accretion rate of 
the SMBHs for each galaxy, we calculate the bolometric luminos¬ 
ity of the AGN as L agn = t]Mbhc 2 where c is the speed of light 
and we assume the radiative efficiency 77 = 0.1. The direct relation 
between LAGN,boi and L([0 II]) is unclear, so we use the results by 
|Heckman et al.| ( [2004] > and |Netzer| ( |2009| > where the [O m]A5007 is 


used as a proxy for the bolometric luminosity of the AGN LAGN,boi 
and the average relation is shown as lAGN,boi/^AGN,[o m] ~ 
3500. Since the line ratio between [O II] and [O III] in AGNs ranges 
from 0.1-5 depending on the types ( [Blandford et al.|l990| , we can 
roughly obtain the fraction of the L([0 II]) out of the bolometric 
luminosity as 3 x 10 -5 < L([0 ll])AGN/LAGN,boi ^ 0.001. 

2.2.6 Dust extinction effect 

In this section, we briefly discuss the effect of the dust extinction 
on the L([0 II]) even though we rather use the intrinsic L([0 II]) 
without dust correction in this study. It is extremely challenging 
to understand the effect of the dust extinction on nebulae emission 
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Figure 6. Fitting of LFs of the [O n] emitting galaxies at z = 0.1, 0.3, 0.6, 1.0, 2.0, and 3.0. For each redshift, we fit the LF with a single power-law at low 
luminosity and Schechter function at high luminosity. We use L([0 n]) = 10 42 0 erg s -1 as a dividing luminosity of the fitting functions except for z = 1.0 
and 2.0 where we apply L([0 n]) = 10 415 erg s -1 , which are shown as dotted vertical lines in each panel. 


lines coupled with the star formation history, metallicity, and the 
evolution as a function of redshift ( [Wilkins et al.|2013b) . 


Both the stellar continuum of the SEDs and emission lines 
are affected by the dust extinction. In observation, the internal dust 
extinction should be corrected to obtain the intrinsic flux at the 
wavelength where the emission lines are located. On the simulation 
side, dust-extinction effect must be added to make a direct compar¬ 
ison with observation. Since the wavelength A = 3730 A is in 
U band, we can apply E(U — V) sta r = 1.0AE{B — V) star to 
the stellar continuum. In galaxies actively forming stars, the em¬ 
pirical relationship of the internal dust extinction ( |Calzetti||2001| ) 
shows that the emission from the ionized gas suffers about twice 
as much reddening as the stellar continuum E(B — V) s t & r = 
0.44^(5 — V)g as . Then the absorption in the U band becomes 
^4u,gas = ^v,gas + (1 MA V ) / (0 A4R V ) . Applying the com¬ 
monly accepted reddening coefficient Rv ~ 3.1 and the extinction 
value at V band Ay — 0.4 mag from which Av, g as = 0.91 is 
obtained [e.g., An a = 0.7^q; 4 for z ~ 0.5 galaxies (Ly et al. 
|2012| >], the extinction coefficient of [O II] emission line is esti¬ 
mated as A [o n] ~ 1.39 mag. Therefore, the dust absorption re¬ 
duces the observed [O II] line luminosity L([0 ll]) 0 to ~ 0.28 
times of the intrinsic flux L([0 ll])i. Many studies find that the 
dust reddening E(B — V) also depends on the intrinsic [O II] lu¬ 
minosity L([0 II] )i. For example, | Kewley, Geller & Jansen| ( [2004) 
measure E(B — V) — (0.174 d= 0.035)log L([0 ll])i — 6.84, so 
the intrinsic [O II] line luminosity is expressed as L([0 II] )i = 
3.11 x 10 -20 L([O II])o ' 495 where the luminosities are given in unit 
of ergs per second. 


3 RESULTS 

3.1 [Oil] luminosity function at low redshift 

We compare the LFs of the [O II] emitting galaxies in the 
MBII simulation with observations at redshift z =0.1, 0.3, and 
1.0 from top to bottom in Fig. [5] Circles show the LFs from the 
MBII simulation while lines show the observational results (Zhu, 


Moustakas & Blanton 2009; Gilba nk et al.|2010] |Ciardullo et al. 


2013[ |Comparat et al.|20l5] > at each redshift. Big circles show all 
the [O II] -selected galaxies while small circles with the Poisson 
errors show the galaxies with AGNs excluded to avoid a possi¬ 
ble AGN contamination. We compare the L([0 II]) with the esti¬ 
mated L([0 II]) agn due to AGN activity for each galaxy using the 
method described in Section [2.2.5 1 We assume that L([0 ii])agn= 
0.001 xLAGN,boi to be conservative. Note that the level of the AGN 
contamination in the L([0 II]) LF is negligible as discussed in Sec¬ 
tion [2]23] 

At z — 0.1, red (observed) and blue (dust-corrected) solid 
lines show HETDEX pilot survey ( [Ciardullo et al.|20L3] ) while the 
red dashed line (observed) shows the result from SDSS ([Gilbank 
|et al.|2010| . The HETDEX pilot survey has a higher number den¬ 
sity compared to SDSS at z = 0.1, but the excess can be explained 
by the cosmic variance caused by the small sample size (< 300 
galaxies at 0.1 < z < 0.56) as discussed in |Ciardullo et al.| ( [2013| ). 
|Comparat et al.| ( [20T5] ) matches well with other observations at low 
L([0 II]) although it fits better with a Schechter function when 
the high L([0 II]) samples are considered together. LFs from the 
MBII are a good matc h with the observation ( Gilbank et al. 2010, 
Comparat et al.|2015|) in the range of [O II] luminosity 10 4u ' b < 
t /in tt1\^ i n 41 - ( j — 1 t in a/trtt cU ^ 7 ™ 


L([0 ii])< 10 41b ergs . The LFs of the MBII show an excess 
in the prediction of bright [O II] emitters at L([0 II]) > 10 41 ' 6 
erg s -1 at low redshift (z ^ 0.3). Note that the galaxy stellar mass 
function (GSMF) of the MBII at low redshift is known to show 


© 2015 RAS, MNRAS 000,[l]{T2] 
























































LF of [Oil] ELGs in the MBII simulation 7 




Figure 7. Evolution of the LF of the [O n] emitting galaxies from z = 3 to 
z = 0.1. In general, the LF increases and the slope of the LF also becomes 
steeper with increasing redshift. There is strong evolution of the LF for 
z > 1.0 but not below z < 1.0. 


an excess compared to observations both in the low- and high-mass 
ends fKhandai et al.|2015) . Since a high fraction of the star-forming 
galaxies in the high-mass end is selected as [O II] emitting galax¬ 
ies, the overabundance in GSMF at the high-mass end is inevitably 
transferred to the [O II] LF. 

At z — 1.0, the LF of the MBII shows a good match with 
|Zhu, Moustakas & Blantonj ( |2009| > throughout the luminosity that 
the DEEP2 survey covers while the LF matches with Comparat 


et al. 


( 2015]) below L([0 ll])=10 43 ° erg s 1 , as shown at the bot¬ 


tom panel of Fig. [5] 


3.2 Fitting function of the [Oil] luminosity functions 

Despite the fact that the [O II] LF is not in perfect agreement with 
observations at low redshift, the [O II] LF at z=l still shows a rea¬ 
sonably good agreement with observations. Therefore we fit the 
[O II] LFs from the MBII for low redshift z ^ 1.0, and extend 
the fitting for LFs at high redshift in the range 1.0 < z ^ 3.0. 
Since the exponential decline of the number density is obvious to¬ 
ward the high luminosity end, we use two fitting functions for lower 
[O II] luminosity than L([0 II]) = 10 42 ergs -1 (dotted vertical 
lines) and high luminosity (L([0 II]) > 10 42 erg s -1 ), respectively. 
For z = 1.0 and 2.0, L([0 II]) = 10 41 ' 5 erg s -1 is used instead as 
the border line to adopt the fact that the Schechter function extends 
to lower luminosity. At low luminosity, single power-law function 
is used while the common fSchechter || 1976| function is used at high 
luminosity. 

The Schechter function has the form 

4>{L)dL = 4>* ex P ( _ ^) d (^) (3) 

where a is the slope of the faint end of the LF, L* is the charac¬ 
teristic luminosity, and 0* is the density of galaxies per magnitude 
with L > L*. At low luminosity (L([0 II]) < 10 42 erg s -1 ) in each 


Figure 8. The slope of the LFs of the [O n] emitting galaxies as a function 
of redshift for the bright (squares) and faint (circles) ends. For both ends, 
a shows a minor change at z ^ 0.3. For the bright end, a increases ap¬ 
proximately as (z + 1) _1 (dot-dashed line) while a of the faint end shows 
the bigger change being proportional to ~ 3(z + l) _1 (dotted line) at high 
redshift z ^ 0.6. The slopes for both ends display the minima at z = 2.0. 


redshift, a power-law fitting function is used with the form 

0(logL)dlogL = 10 ( “ i+1)(losZ '“ 42 ' 0)+/3i dlogL (4) 

while the Schechter function is used for high luminosity as 

0(logL)dlogL — io ( cx h+ 1 ) l °g(ir*)+f 3 h eX p (— ) dlogL. (5) 

Jj* 

Fig. [6] shows the LFs along with fitting parameters at redshift 
z — 0.1, 0.3, 0.6,1.0, 2.0, and 3.0, respectively. We apply Poisson 
errors to estimate the uncertainty of the number of galaxies in each 
luminosity bin (see Appendix[B]for the comparison of the estimated 
Poisson and Jackknife errors). Table[l]shows the fitting parameters 
and the corresponding reduced % 2 values at different redshifts. At 
high luminosity ends, the number of galaxies are dominated by the 
Poisson errors. For example, only 2-3 galaxies are found in the sim¬ 
ulated volume at the luminosity bin of L([0 II]) ~ 10 43 6 erg s -1 . 
However, at the low luminosity ends, the small Poisson errors due 
to the high number of galaxies (e.g., ~ 73000 galaxies at the lumi¬ 
nosity bin L([0 II]) ~ IO 40 ' 8 erg s -1 for z — 3.0) return the high 
values of x 2 - 

Fig.[7]shows all the [O II] LFs from z — 3.0 to the local uni¬ 
verse z = 0.1 to show the LF evolution as a function of redshift. In 
general, the number of galaxies increases and the slope of the LFs 
decreases (i.e., a decreases) as a function of redshift. The number 
of [O II] galaxies at the lowest luminosity increases by almost 2 
orders of magnitude from z = 0.1 to z = 3.0 whereas relatively 
small change is observed at the highest luminosity end. The LFs 
do not evolve much at low redshift z < 1.0 but strong evolution is 
observed at high redshift 1.0 ^ z ^ 3.0. 

Fig. [8] shows the slope a of the LFs as a function of redshift. 
The slopes for the bright (L([0 II]) ^ IO 42 0 ergs -1 ) and faint 
(L([0 II]) ^ 10 41 ' 2 ergs 1 ) ends of the LFs are plotted against 
(z + l) -1 as squares and circles, respectively. For both ends, a 
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Table 1. Fitting parameters of the LFs of [Oil] emitters in equations j4j and |5j 


z 

on 

A 

W,red 

Oih 

p h 

L* 

^-h,red 

0.1 

-2.21 ± 0.06 

46.96 ± 2.60 

1.10 

-1.59 ±0.15 

-4.10 ±0.22 

43.06 ±0.14 

0.16 

0.3 

-2.23 ± 0.07 

47.88 ± 2.90 

2.16 

-1.91 ±0.13 

-4.80 ± 0.34 

43.40 ± 0.21 

0.14 

0.6 

-2.10 ±0.08 

42.64 ± 3.22 

3.37 

-1.88 ±0.14 

-4.56 ± 0.30 

43.16 ±0.18 

0.09 

1.0 

-2.40 ± 0.26 

55.37 ± 10.77 

38.03 

-2.05 ±0.11 

-4.13 ±0.34 

42.98 ± 0.21 

1.49 

2.0 

-3.15 ±0.14 

86.61 ± 5.92 

37.80 

-1.74 ±0.06 

-3.35 ±0.14 

42.85 ± 0.09 

1.14 

3.0 

-2.59 ± 0.07 

64.10 ±3.07 

43.42 

-1.65 ±0.08 

-3.29 ±0.11 

42.96 ± 0.07 

0.19 



log(L[OIII] erg s 1 ) 


Figure 9. LFs together with fitting functions of the [O ill] ELGs at z = 0.1, 0.3, 0.6, 1.0, 2.0, and 3.0. The [O ill] LF is similar to the [O n] LF except for 
the fact that L([0 ill]) is higher than L([0 n]) of a same SFR. The dotted vertical line at each redshift shows the dividing luminosity for fitting LFs at low and 
high luminosity. 


shows a minor change at z ^ 0.3. For the bright end, the slope 
changes approximately as (z + l) -1 whereas the slope of the faint 
end shows the bigger change approximately being proportional to 
~ 3(z ± l) -1 at high redshift 0.6 ^ z ^ 2.0. The slopes for both 
ends display the minima at z = 2.0. Since the Poisson errors dom¬ 
inate at the bright end of the LFs, the slopes for L([0 II]) ^ 10 42 ' 4 
ergs -1 are plotted additionally (pentagons) to see how the errors 
affect the determination of the slopes. The number of [O II] galax¬ 
ies at the highest luminosity bins L([0 II]) ~ 10 43 ' 6 ergs -1 is 
typically 2-3 in the simulated volume, and thus the Poisson statis¬ 
tics greatly affects the determination of the slopes. For a compari¬ 
son with observations, the slopes from|Zhu, Moustakas & Blanton| 
( |20Q9| >; [Gilbank et al. 1(2010) ; |Comparat et al.| ( |2015| ) are plotted 
together in Fig. [8] however, note that their luminosity cuts do not 
exactly match with the ones from the current work for both bright 
and faint en ds. F or e xampl e, the luminosity range of the bright end 
from |Gilbank et al. 1(2010) is 10 41 0 < L([0 II]) < 10 42 5 ergs -1 . 
The evolution of the [O II] LF as a function of redshift is driven by 
the evolution of GSMF, i.e., by the cosmic SFR in different stellar 
masses. Thus, it is important to model the SFR in various stellar 
mass ranges. At low mass end, the SFR is greatly affected by the 


star-formation feedback model whereas the AGN feedback greatly 
affects the high stellar mass end. 


3.3 [OIII] luminosity functions 

In this section, we briefly discuss the LFs of the [O III] emit¬ 
ting galaxies in the MBII, since the upcoming survey WFIRST- 
AFTA will explore the evolution of the [O III] emission line LFs at 
high redshift. We repeat the process for the LFs of the [O II] emit¬ 
ters for [O III] emission line. In the synthesised SEDs of the 
star-forming galaxies, the emission-line fluxes of the [O II] and 
[O III] lines are calculated in a similar manner. Note that the 
[O III] LFs in Fig. [9] are very similar to the [O II] LFs in Fig. [6] 
except for the fact that the [O III] luminosity is higher than [O II]. 
We fit the LFs with two functions as we do for [O II] LFs; how¬ 
ever, we split the LFs at the luminosity 0.4 dex higher than the 
luminosity for the [O II] at L([0 III]) = 10 42 ' 4 ergs -1 . Fig. [ 9 ] 
shows the [O III] LFs along with fitting functions at redshift z — 
0.1,0.3,0.6,1.0, 2.0, and 3.0, respectively. Fig. [l()| shows all the 
[O III] LFs at 0.06 ^ z < 3.0 together; Table |2jlists the fitting 
parameters. 
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Table 2. Fitting parameters of the LFs of [OIII] emitters in equations j4j and |5j 


z 

ai 

A 

^l,red 


p h 

L* 

^•h,red 

0.1 

-2.23 ± 0.08 

48.10 ±3.51 

0.95 

-1.90 ±0.26 

-4.43 ± 0.55 

43.50 ± 0.33 

0.24 

0.3 

-2.22 db 0.06 

47.88 ±2.51 

0.69 

-2.24 ±0.16 

-5.76 ± 1.13 

44.22 ± 0.71 

0.17 

0.6 

-2.01 ± 0.08 

39.57 ± 3.52 

1.75 

-2.01 ± 0.20 

-4.55 ± 0.44 

43.43 ± 0.24 

0.09 

1.0 

-1.80 ±0.04 

30.96 ± 1.49 

0.20 

-2.22 ± 0.08 

-4.09 ± 0.25 

43.23 ±0.13 

0.20 

2.0 

-2.80 ±0.18 

73.08 ± 7.36 

11.52 

-1.67 ±0.10 

-3.24 ±0.16 

43.16 ±0.10 

1.28 

3.0 

-2.86 ± 0.04 

75.79 ± 1.52 

3.08 

-1.67 ±0.07 

-3.23 ± 0.08 

43.23 ± 0.05 

0.10 



Figure 10. Evolution of LF of the [O in] emitting galaxies from z = 3 
to z = 0.1. The overall shape of the [O ill] LF is similar to the one for 
[O II] LF with a shift in luminosity. 

3.4 Evolution of the auto-correlation function of 
[Oil] emitters 

In this section, using the spatial distribution of the [O II] emission¬ 
line selected galaxies, we explore the evolution of the 2-point 
auto-correlation function [£(r) =DD(r)/RR(r) -1] at 0.06 ^ z ^ 
4.0, where DD(r) is the number of galaxy pairs with separa¬ 
tion r and RR(r) is the number of pairs with the same sepa¬ 
ration for a random (i.e., Poisson) distribution. Fig. El shows 
the auto-correlation functions of [O II] emitters at the redshifts 
z — 0.06,0.1,0.3,0.6,1.0, 2.0,3.0, and 4.0 along with the cor¬ 
responding fitting functions shown as dashed lines. When the 
auto-correlation function is fitted with a power-law as £fit(r) = 
(r’// C orr) mcorr , £fit(fcorr) = 1 is obtained where l co rr is the cor¬ 
relation length within which the galaxy distribution is correlated. 
In Fig. m the intersection between the power-law fitting func¬ 
tion for each redshift £fi t (r) and £ = 1 shows the correlation 
length l corr for each redshift. Table [3] lists the fitting parameters 
of the auto-correlation functions. The £ shows a significant evolu¬ 
tion from z = 2.0 to 1.0 whereas it evolves mildly at z < 1.0 
or z ^ 2.0. At z < 1.0, the slopes remain constant in the range 
m co rr = —1.55 to —1.45 while the correlation length increases 
from l corr = 4.47 Mpc h~ x at z = 1.0 to / corr = 5.23 Mpc h~ x 
at 2 ? = 0.06. At z ^ 2.0, the slopes display slightly larger values 



log(r [Mpc h 1 ]) 


Figure 11. The auto-correlation function [£(r) =DD(r)/RR(r) -1] of the 
[O II] emitting galaxies for redshift z = 0.06, 0.1, 0.3, 0.6, 1.0, 2.0, 3.0, 
and 4.0. £ shows a significant evolution from z = 2.0 to 1.0. The inter¬ 
section between the power-law fitting function £fi t (r) (dashed) and £ = 1 
(dotted) shows the correlation length l COT r for each redshift. 


Table 3. Fitting parameters of the auto-correlation functions 


z 

Ffl corr 

lcorr (MpC h 1) 

0.06 

-1.52 

5.23 

0.1 

-1.55 

4.99 

0.3 

-1.45 

4.66 

0.6 

-1.52 

4.92 

1.0 

-1.50 

4.47 

2.0 

-1.38 

3.37 

3.0 

-1.33 

3.10 

4.0 

-1.35 

3.31 


than the ones for z ^ 1.0 in the range of m co rr = —1.38 to —1.33 
and the smaller correlation lengths are clearly observed in the range 
lcorr = 3.10-3.37 Mpc hr 1 . 


4 SUMMARY AND DISCUSSION 

In this paper, we investigate the properties of the [O II] ELGs 
in the state-of-the-art high-resolution cosmological simulation 
MBII. The MBII simulates galaxies including baryonic physics 
for star-formation and AGN feedback in a comoving volume of 
100 h~ x Mpc on a side from z — 159 to z = 0.06. From the 
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synthesised SEDs of the individual galaxies, which includes the 
stellar continuum and the emission lines based on the star forma¬ 
tion history, we select a sub-sample of star-forming galaxies based 
on the [O II] luminosity in the redshift range 0.06 ^ z ^ 3.0. We 
validate the use of the extracted L([0 II]) by comparing it with sev¬ 
eral observations. The MBII simulation shows a good agreement 
with observations, and we focus on the evolution of the [O II] (and 
[O III]) LFs in the redshift range 0.06 < z ^ 3.0. We expect that 
the current work is useful for forthcoming surveys such as HET- 
DEX, DESI, and WFIRST-AFTA. We summarise our main findings 
as follows. 


• The specific SFR as a function of stellar mass agrees with 
the GAMA survey when galaxies with high [O II] luminosity are 
selected as L([0 II]) ^ 10 41 ' 5 erg s -1 . 


• We show that the [O II] LF at z ®= 1.0 from the MBII shows 
a good agreement with the LFs from several surveys ( |Zhu,| 
|Moustakas & Blanton| |2009| |Comparat et al.| |2015| > below 
L([0 ll])=10 43U erg s -1 while the low redshifts (z < 0.3) show an 
excess in the prediction of bright [O II] galaxies, but still displaying 
a good match with observations ( |Gilbank et al.| |2010, Ciardullo 
et al.|2013~ Comparat et al.|2015| below L([0 ll])=10 4i ^ ergs -1 . 


• We present the LFs of [O II] and [O III] ELGs at the redshift 
range 0.1 ^ z ^ 3.0 and provide fitting functions for each redshift. 
Each [O II] LF at different redshift is fitted with a single power-law 
at low luminosity (L([0 II]) < 10 42 0 ergs -1 ) while Schechter 
function is applied at high luminosity (L([0 II]) > 10 42 0 erg s -1 ). 
In general, the LF increases and the slope of the LF also becomes 
steeper with increasing redshift. The slopes of the LFs at bright 
and faint ends range from -3 to -2 showing minima at z = 2. The 
slope of the bright end evolves approximately as (z + l) -1 at 
z ^ 2 while the faint end evolves as ~ 3 (z + l) -1 at 0.6 ^ z ^ 2. 
We apply a similar analysis to the [O III] LF with a shift of 0.4 dex. 


• The auto-correlation function of [Oil] ELGs shows a signifi¬ 
cant evolution from z = 2 to 1 while it changes mildly at z ^ 1 or 
2.0 ^ z ^ 4.0. The correlation length increases from ~3 Mpc /i -1 
for z ^ 2 to ~5 Mpc h~ x for z ^ 1. 


The current theoretical study of the LFs of [O II] (and [O III]) 
emitting galaxies will be useful for the forthcoming surveys HET- 
DEX, DESI, and WFIRST-AFTA. The current high resolution of 
MBII simulations makes it possible to investigate the [O II] line 
of an individual galaxy and construct a better theoretical model for 
LFs, but the low luminosity end still suffers from shot noise. We ex¬ 
pect that future simulations with even higher resolution can provide 
a solid feature at low luminosity. At the bright end of the LFs where 
Poisson statistics dominates, enlarging the simulation volume will 
help to get better errors for each luminosity bin. 

The excess of galaxies at high luminosity end in [O II] (and 
[O III]) LF is inevitably transferred from the overabundance at the 
high mass end in GSMF of the MBII. In this respect, the future 
theoretical study of [O II] LF should be focused on the GSMF by 
applying observation matching star-formation and AGN feedback 
models. Especially at the faint end, a stronger star-formation feed¬ 
back model (e.g., Okamoto et al.|2010| used in Illustr is ( [Vogels- 1 
|berger et al.|2014| and BlueTides ( [Feng et al.|2015| ) can reduce 
the number of star-forming galaxies. On the bright end, increasing 
the AGN feedback (e.g., EAGLE simulation by |Schaye et al.|2015) 
can suppress the star-formation by driving the gas out of the halos. 


The dust obscuration in star-forming galaxies in low and high 
redshift galaxies remains as an open question. Future studies should 
consider the effect of dust reddening on the LF at different red- 
shifts. 
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APPENDIX A: VALIDATION OF THE [Oil] LF 

We compare the [O II] LFs calculated using different methods for 
the purpose of validation in Fig. |A1| We confirm that our method 
recovers the intrinsic LF of the [O II] well, as shown in Fig |Al| The 
sampling of the [O II] emitting galaxies start to lose star-forming 
galaxies at L([0 II]) ^ 10 4O ' 6 ergs -1 due to the sensitivity of the 
[O II] emission line on the SF timescale and the shot noise due to 
the mass resolution of the gas particles in the MBII simulation. 



log(L([Ollj) erg s 1 ) 


Figure Al. LFs of [O II] emitting galaxies in the MBII at z = 0.06. Cir¬ 
cles show the [O n]-selected galaxies while squares show the LF of the 
raw L([0 n]) before the emission lines are added to the stellar SEDs for 
each galaxy. Triangles show the LF of the L([0 n]) which is converted 
from the SFR using the empirical relationship by [Kewley , Geller & Jansen| 
([20041. Our method of selecting [O II] ELGs recovers the LF of the intrin¬ 
sic [O II] lines well. The LF starts to lose star-forming galaxies due to the 
fact that the [O II] lines trace the average SF history of the past ~ 20 Myr. 
Note that the LF converted from the SFRs keeps increasing toward the low 
luminosity end. Shot noise due to the mass of the star particles also starts 
to dominate in the low luminosity L([0 n]) < 10 40-6 ergs _1 shown as a 
vertical dashed line. 
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Figure A2. Comparison of Poisson and Jackknife re-sampling errors. For 
each luminosity bin, Poisson errors (VN) are obtained (circles) from the 
number of galaxies ( N ). Jackknife sampling of 8 cubes with comoving vol¬ 
ume of 50 Mpc h~ x on a side is used to get errors due to the cosmic vari¬ 
ance. The Poisson and Jackknife errors, shown here, are normalised by the 
mean value. 
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APPENDIX B: STANDARD ERRORS OF THE LF 

We compare the errors of the LFs in each redshift from Poisson 
statistics and Jackknife re-sampling. Fig. |A2| shows the standard 
errors normalised by the mean values for each luminosity bin. Pois¬ 
son errors (y/N) are obtained (circles) from the number of galaxies 
(N) in each luminosity bin. Jackknife sampling of 8 cubes with 
comoving volume of (50 Mpc /i -1 ) 3 is used to get the errors to 
consider the cosmic variance. Poisson errors are bigger than the 
Jackknife re-sampling errors, whereas both errors are comparable 
at low luminosity L([0 II]) ~ 10 41 erg s -1 . 
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